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ABSTRACT 

GD-1 is a long, thin, Milky Way star stream that has readily visible density variations along its 
length. We quantify the locations, sizes and statistical significance of the density structure, i.e. gaps, 
using a set of scaled filters. The shapes of the filters are based on the gaps that develop in simulations 
of dark matter sub-halos crossing a star stream. The high Galactic latitude 8.4 kpc long segment of 
GD-1 that we examine has 8 ± 3 gaps of 99% significance or greater, with the error estimated on the 
basis of tests of the gap-filtering technique. The cumulative distribution of gaps more than 3 times the 
width of the stream is in good agreement with predictions for dark matter sub-halo encounters with 
cold star streams. The number of gaps narrower than 3 times the width of the GD-1 stream falls well 
below the cold stream prediction which is taken into account for the gap creation rate integrated over 
all sizes. Simple warm stream simulations scaled to GD-1 show that the falloff in gaps is expected for 
sub-halos below a mass of 10 6 M Q . The GD-1 gaps requires 100 sub-halos > 10 6 M within 30 kpc, the 
apocenter of GD-1 orbit. These results are consistent with LCDM sub-halo predictions but further 
improvements in stream signal-to-noise and gap modeling will be welcome. 
Subject headings: dark matter; Local Group; galaxies: dwarf 



1. INTRODUCTION 

The GD-1 star stream was discovered in the Sloan 
Digi tal Sky Survey (SPSS) Da ta Release 4 photome- 
try (| Grillmair fc Dionatosl [2006) with Data Release 7 
improving both the sky coverage and photometr ic uni- 
formity ^metEeEaD iSffil EoEoiovSan[20i3). The 
stream arcs more than 80° across the northern sky, pass- 
ing within 30° of the Galactic pole. The visible section of 
the stream is about 8 kpc distant at its midpoint with the 
distance increasing from about 7 to 11 kpc. The stream 
is exceptionally narrow, having an angular width of 0.5°, 
which corresponds to a linear width of about 70 pc. The 
stream has a remarkable length to width ratio of at least 
100. Figure [1] rotates the sky to place the stream near 
the reoriented equ ator in a Mercator projection. 

SEGUE spectra (jYannv et al. 2009) provide velocities 
and metallicity measurements of the stars in selected re- 
gions along the s tream. Those velocities, augmented with 
proper motions (|Munn et al.l 120041 120081) f rom the com- 
parison of the SPSS to the USNO catalogs (jWillett et al.l 
l2009t IKoposov et al.ll2010l ). provide sufficiently accurate 
phase space information along the stream to give a good 
GD-1 orbit and some limits on the shape of the Galac- 
tic potential. The derived orbits have a perigalacticon 
of about 14 kpc and apogalacticons of 26-29 kpc. The 
orbits passes through the plane of the galaxy at large 
radii where there are relatively few high density molec- 
ular gas clouds or HII regions. The narrowness of the 
stream and the stella r mass in the stream, 2 x 10 4 M Q 
(jKoposov et alJ I2010I ). suggest th at the likely progeni- 
tor is (or was) a globular cluster ([Grillmair fc Dionatosl 
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2006). No progenitor for the GD-1 stream has yet been 
identified close to the orbit of the stream. 

The stream has readily visible density variations along 
its length, whi c h are the primary interest of this paper. 
IKoposov et al.l (|2010D simply stated that the reason for 
the density variations is not clear. For our purposes, 
an absent or distant progenitor is a useful simplifica- 
tion in that the epicyclic d ensity variations discussed in 
iKupper et aT] (|200l. I2012T) play no role in the density 
variations of the stream. 

The GD-1 stream allows a fairly strong test for the exis- 
tence and properties of the large population of sub-halos, 
N(> M) ~ 1.6 x 10 5 (M/10 6 M Q )-°- 9 , orbiting within 
a Milky Way-like dark halo (counted to a distance of 
433 kpc) predicted to exist in dark halos simulated from 
LCDM initial conditions (iDiemand. Kuhlen fc Madaul 
[20071: ISprmgel et all 120081 IStadel et al.ll2009fl . When a 
sub-halo passes through a stellar stream it normally pulls 
out a loop of sta rs which develops into a visible gap in the 
strea m density (fCarlbcrg 2 0091 : lYoon. Johnston fc HoggJ 
120111 ). The rate per unit length of gap creation, R u , is de- 
rived in the col d stream approximation as i?u M~ 36 
(Carlberg 2012), where M m is the smallest mass sub-halo 
that can create a visible gap of any size greater than the 
stream width. 

We previously developed a density filtering procedure 
to locate gaps and assign a statistical confidence level 
to their detection wh ich we applied to the Pal 5 stream 
([Carlberg et al.l I2012T ). Now that improved GD-1 data 
are available it is a better case for the study of gaps. 
That is, Pal 5 is significantly more distant than GD-1, 
some 23 kpc as compared to 8 kpc, respectively, and, 
the Pal 5 stream skirts the Galactic plane in projection, 
hence has a fairly large stellar background. Furthermore, 




Fig. 1. — An image of the stream in a rotated sky coordinate system in which the stream is placed near the new equator. The coordinates 
are in degrees relative to the second cross-hair on the figure and are not aligned with any standard system. To remove the varying 
background, the masked image has been smoothed over 5°, subtracted from the original image, and then smoothed with a 0.4° Gaussian. 
The "cross-hairs" are used as the initial guess of the centerline used to extract the stream. 



the Pal 5 globular cluster is very much present with the 
galactic tides still pulling stars away from the cluster 
on orbits that are expected to crea te density variations 
in the stream near to the cluster (jDehnen et al.l 120041 : 
IKupper et al.l 120121) . The GD-1 star stream is a high 
galactic latitude system, nearly a segment of a great cir- 
cle when viewed in re-oriented coordinate s (see Figure [1]) 
(jWillett et all 120091: IKoposov et""aHl2010[ ). It has a rela- 
tively high signal-to-noise along its length and is one of 
the best streams to search for gaps and to test theories 
for their evolution in the galactic potential. 

In this paper, we first search for the gaps along the GD- 
1 stream using a set of filters based on the density profiles 
of gaps that develop in simulations of sub-halo stream 
interactions. We characterize the statistical behavior of 
the filters on a model stream with gaps inserted. The 
numbers of gaps in the GD-1 stream and the distribution 
of their sizes are compared the predictions from dark 
matter sub-halos to estimate the lowest effective sub-halo 
mass and their numbers within the orbit of the stream. 

2. THE STREAM MAP AND DENSITY PROFILE 

The SDSS Data Release 8 (lYork et al.l \2(M 
lAihara et aLll2011h extinction corrected stellar photome- 
try is optimally match-filtered in color-magnitude space 
to identify old metal poor stars at distances around 
8 kpc. The filte r ing u se s the procedures documente d 
in iRockosi et al.l (120021); iGrillmair fc Dionatosl (|2006l) : 
iGrillmairl (|2009t ) and lGrillmairl (|2011[ ). The pixel values 
are star counts filtered on the basis of their agreement 
with the M13 color-magnitude relations and weighted 
with a luminosity function for a Z = 0.0003 low metal- 
licity population. 

The distance to the GD-1 stream varies along its 
length. The same color-magnitude filter has been offset 
+0.2 magnitudes to make four separate sky maps with 
approximate distances ranging from 6.7 kpc for the first 
plane to 8.8 kpc for the last (based on an assumed dis- 
tance of 7.7 kpc for M13, the basis of the filter). We 
will present the density profiles of all four maps below. 
Most of the results below use the 8.1 kpc distance version 
(the third of the four maps) on the basis that it has the 
highest total stream density. 

We rotate the map from its J2000 Mercator coordinates 
to a new Mercator projection with the stream along the 
new equator. The (arbitrary chosen) stream point orig- 
inally at [RA, Dec] = [141.80°, 27.04°], where the sec- 



ond cross-hair symbol is located, becomes the new [0,0] 
and the new equator passes through the point [RA, Dec] 
= [222.29°, 58.22°]. The rotated image is placed into the 
same size pixels, 0.1°, as the original image, which in- 
troduces correlations between adjacent pixels, but the 
pixels remain statistically independent at larger separa- 
tions. Figure [1] shows the image in the new coordinate 
frame. For better viewing a 5° Gaussian smoothed im- 
age is subtracted to approximately remove the variation 
in the stellar density over the sky. The circled plus signs 
show the points that serve as our initial estimates of the 
centerline of the stream. Although the path of the stream 
can be identified further after it crosses the Sagittarius 
stream (the broad feature that runs through [-5,0] on the 
left side of the map) we do not include that portion in 
our analysis since it has undesirably large errors from the 
large and variable background correction. On the right 
side it is also unclear exactly where the stream goes, al- 
though we do note that our 70° point extends the stream 
a few degrees further than previous papers. 

The segment of the GD-1 stream that we have selected 
is ab ove 40° galact i c latit ude, where the E(B-V) values 
from lSchlegel et all ()1998h are generally below 0.02 mag. 
Any extinction corrections that equally affect the stream 
and background region would make little difference to 
the background subtracted stream density. Because the 
total extinction corrections are so low, 5% or less, in a 
stream of about 20% over-density, any errors that may 
exist in the extinction corrections are not sufficient to 
cause gap-like density variations. 

The stream is given xy coordinates relative to the cen- 
terline. The initial placement of a set of centerline points 
is done by hand on local density peaks, after which a 
simple random search is undertaken to maximize the 
mean stream density over the interval [10,70]. The op- 
timization procedure moves the points less than a pixel 
on average and boosts the density relative to the ini- 
tial guess by about 10%. A spline interpolates the cen- 
terline between the points. The x coordinate is mea- 
sured as distance along the fitted centerline, relative to 
the first point which is at the edge of the Sagittarius 
stream. The y coordinate is the distance to the interpo- 
lating spline. The centerline is a small amplitude, slowly 
varying curve. However there are some small scale local 
excursions which are not removed. 

Perhaps the most interesting local excursion is located 
near x = 50° in Figure [2j where the stream has an s-like 
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Fig. 2. — An image of the extracted stream relative to the optimized centerline. The image is smoothed with a 0.2° Gaussian. The 
coordinates are in degrees relative to the centerline with x and y and coordinates being along and transverse to the stream, respectively. 
No compensation for the varying background has been applied. Our analysis is applied to the [10,70] degree range of this image. 
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Fig. 3. — The density profile across the stream, averaged from x 
= 10-70°. The FWHM is 0.50°. The stream is extracted within 
±0.25° of the fitted centerline and the background region starts at 
±0.5°. 



shape above and below the axis. The feature is visible 
in Figure [1] around 42°of its rotated coordinate system 
(which is not the same as the Figure [5] system) . The 
feature, although not very compelling at the available 
signal-to-noise, could be the result of an encounter with 
a sub-halo or possibly be the location of a recently dis- 
rupted progenitor object. We consider an alternate cen- 
terline through these points below. 

We estimate the background using the regions on both 
sides of the stream beyond ±0.5° where no stream signal 
is visible. The background region extends 5 times the 
width of the selection region, that is, from 0.5 to 1.75° 
distant from the centerline on either side of the stream 
in the image with no smoothed background subtracted. 
This background is subtracted from the total density in 
the centerline region to yield the linear density of the 
stream. The background subtracted density is extracted 
into bins of 0.1° along the stream which we use in the 
analysis of the GD-1 stream properties. 

The stream densities are summed along its length from 
10-70° in Figure [5] to create a transverse density profile. 
We measure the FWHM of the stream (simply using the 
full width at half the peak height) over the 60 degree 




Fig. 4. — The background subtracted density along the stream 
in 1° bins with distance increasing from the bottom upwards, sub- 
sequent lines offset by 0.1 in density. The dashed line indicates the 
estimated error in each bin. 



length to be 0.50°, as shown in Figure [3] Koposov ct al. 
(2010) over a similar length of the stream found a Gaus- 
sian width of 0.20°, which, for a Gaussian, converts to a 
FWHM of 0.47°, which is consistent with our measure- 
ment. To create a density profile along the length of the 
stream we sum the density over ± 0.25° around the cen- 
terline which captures about 90% of the stream stars. 
Larger stream widths slightly increase the total density 
of the stream, but it leads to a decrease of the signal to 
noise. 

An estimate of the point by point error of the mea- 
sured stream density is calculated from the dispersion in 
the estimate of the local background density, cTj. Since 
the background is N times the width of stream region 
the background noise in the stream region will be ^/N &{ . 
The mean density in the stream is a factor 5 higher than 
the background which boosts the error in the stream 
to <Jiy/N{l + S). Subtracting the background from the 
stream adds its error of Oi in quadrature to the stream er- 
ror to give an error for the background subtracted stream 
to be <r iy /N(l + 8) + 1. The GD-1 stream has 5=0.178 
in our 8.1 kpc map and normally we use a background 
width such that N — 5 although we do the analysis for 
N = 7 as well to check for anomalies. Including the 
stream over-density in the variance calculation boosts 
the required peak height for 99% confidence by about 
5%. 

The density profiles for all four distances of the 
matched filter are displayed in Figure |4j Each has had 
its own search for an optimal centerline and each uses its 
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own background calculation. The profiles are highly cor- 
related. The mean density in the streams is, from nearest 
to farthest, 0.147, 0.170, 0.178, and 0.166, respectively. 
Given the relatively minor differences between the pro- 
files we prefer to use the profile with the largest mean 
density, the one at a mean distance of 8kpc (the second 
from the top in Fig. [3} as the best to search for gaps. 
However, we do analyze all four distance maps below. 

3. STREAM GAP-FINDING 

The passage of a dark matter sub-halo, or in principle 
any massive object, through or close to the stream causes 
material in the region of encounter to gain or lose angular 
momentum and spread apart creating a gap in the mid - 
dle and density pileups on either side (|Carlbergll2012h . 
To find such gaps we use th e density profile filteri ng ap- 
proach initially developed in lCarlberg et al.l (|2012t ). aug- 
mented below with a study of its statistical performance 
for the GD-1 stream parameters. 

3.1. Scaled Gap- Filters 

The first step in gap-finding is to convolve the stream 
density profile with a gap filter which approximates 
the shape of the gaps found in simulations. We use 
gap-filters with the functional forms wi(x) = (a; 6 — 
l)exp (-1.2321a; 2 ) and w 2 (x) = (x 8 - 1) exp (-0.559a; 4 ) 
both of which have a low density floor and compensat- 
ing peaks on either side so that the filters have zero mean 
over the range x = [—3, 3]. The u>2 filter has a quicker de- 
scent to the floor and sharper compensating peaks than 
w\. The filter is scaled to gap lengths of 0.1° to 5° along 
the stream. 

The second step in gap-finding is to establish the con- 
fidence levels of the peaks in the filtered density profile, 
that is, the probability that a peak in the filtered density 
field is a gap detection and is not a peak that results from 
the noise in the stream density measurement. A confi- 
dence level of, say, 99%, means that there is only a 1% 
chance that the peak is the result of chance fluctuations 
in a noisy but otherwise gap-less density profile. 

We use the measured background density profile as 
the starting point to construct random density profiles 
which do not contain any gaps. The background den- 
sity field is first flattened with a third order polynomial 
to remove any gradients along the stream density pro- 
file which can artificially increase the scatter around the 
mean stream density. The flattened background density 
field is then scaled to have the same estimated noise level 
as the stream. This gap-free density field is then sam- 
pled, with replacement, to create 100,000 random den- 
sity profiles which are then analyzed for peaks in pre- 
cisely the same way as the stream. It is this step that 
requires that the background density profile have an ap- 
proximately constant mean density, otherwise in the re- 
sampling points from regions of different mean densities 
are mixed which leads to an increased number of peaks 
being found. The gaps found in the random realizations 
are sorted to determine the filtered field heights required 
to be above some required confidence level. We normally 
use the 99% level in this paper. After applying filters 
of all lengths to the stream density profile we select the 
most significant peak above 67% probability present at 
any location in the stream. This procedure means that 



lower significance narrow gaps are appropriately incorpo- 
rated into higher significance wider gaps. The edges of 
the highest significance peak detection at each location, 
as defined by the filter scaling, is used for the purpose 
of illustration to create a square shouldered gap with a 
depth of 0.5 divided by its confidence level relative to the 
required level as shown in Figures [5] and [7] 

3.2. Gap Filter Tests 

We undertake a number of simple tests of the statisti- 
cal reliability of the gap-finding procedure for a stream 
with the noise characteristics of GD-1. Non-overlapping 
gaps are inserted into a uniform density using our fil- 
ters (normally w\ is used) as shown in the lower line in 
Figure [5] to create an artificial stream with gaps. Gaus- 
sian noise is added at each point as shown as the noisy 
density profile in Figure [5] The measured RMS fluctu- 
ations in the data are 34% of the mean level, but that 
also includes the extra fluctuations from whatever gaps 
are present. The artificial stream has noise added with 
0" = 0.3 around a mean of one. We also create a realiza- 
tion of the background with the same noise level. These 
two profiles are then analyzed in exactly the same way as 
real stream data. The resulting gap detections are shown 
in Figure [5] 

The statistical differences between the distribution of 
input gap sizes and the gap-filter results is quantified in 
Figure |5] which shows the input distribution of gap sizes 
(triangles) and the recovered 95 and 99% confidence dis- 
tributions for the w\ filter. The W2 filter is not well 
matched to the inserted w% gaps and always has more 
false positives than the w\ filter so we do not consider it 
further. The w\ filter performs well for 99% confidence 
gaps even at a noise level enhanced above what we mea- 
sure in the data. The gap distributions from the w± filter 
at a = 0.3 noise for 100 realizations of the added noise 
are shown in Figure |51 On average the total number 
of recovered gaps at 99% confidence level is 1.44 ± 0.15 
times the input number. We will take the mean bias to 
be a factor of 1.44 with a 3 standard deviation (~99% 
confidence) spread of 30%. An extremely important out- 
come is that at small angles the recovered distribution 
is shifted upwards but otherwise retains the same basic 
N(> t) shape with no induced features. 

The biggest gaps are biased to small sizes, but less 
than Figure |6] suggests. That is, the procedure finds 
that the largest recovered gap is on the average 2.4 ± 0.4 
whereas the input value is 3. Consequently the recovered 
largest gaps is on the average at 1.5 standard deviations 
low, but well within the population distribution. Fur- 
thermore we caution that this result is for only a single 
realization of a gap distribution. Further efforts to im- 
prove the performance of the gap-filtering procedure are 
clearly warranted. 

These simple tests also underscore that the quality of 
the results will depend on the shapes of the gaps and 
the details of the noise distribution. We note that our 
assumed gap density profile and the Gaussian noise that 
we have added are only approximate representations of 
the properties of the real data. 

3.3. Gaps in the GD-1 Stream 

The results of gap- filtering the GD-1 stream with the 
w\ filter are displayed in Figure [7] The w\ filter finds 
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Fig. 5. — A demonstration of the gap filtering procedure with 
test data. The noiseless stream is shown at the bottom and the 
same stream with added noise is above it. The 99% confidence 
gaps that the filtering procedure finds are the square profiles that 
are below a value of 0.5 in the plot. 
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Fig. 6. — Comparison of the distribution of gap sizes input (trian- 
gles), and the mean of 100 random distributions at 90 (dotted) 95 
(solid line) and 99% confidence (line with small diamonds) found 
for a = 0.3 with the w\ filter. 

that there are 12 gaps above 99% confidence. Our tests 
of gap-filtering indicates that we should reduce the 99% 
confidence results by a factor of 1.44, to give 8.3 gaps 
and then assign a 30% error to the result, so 8 ± 3 at 
99% confidence. Although we prefer our 99% values it 
is interesting to note that at 95% confidence the same 
procedures find 18 gaps in the GD-1 stream segment. 
The statistical tests find that this is a factor of 2.2 above 
the true number, which then becomes 8, the same as 
estimated at 99% confidence. 

We undertake precisely the same gap analysis on each 
distance slice to check for anomalies and significant vari- 
ations. From nearest to farthest there are 7, 6, 12 and 
11 99% gaps, respectively, with the third distance being 
our preferred choice. At 95% confidence the numbers are 
21, 20, 19, and 22 gaps, respectively, although our test- 



FlG. 7. — The GD-1 reconstructed gaps at 99% confidence with 
distance increasing from the bottom upwards. The density profiles 
in 0.5° bins are also shown. Our primary results are for the 8.1 
kpc distance, second from the top. 

ing indicates that about half of these are false positives. 
However, these results indicate that the streams at each 
distance are strongly correlated and as a consequence, so 
are the gaps that are found. 

The region around 50° of Figure [5] contains an s-like 
feature around 50° which is not traced by our slowly 
varying streamline function. If we force the streamline 
to go through these points it increases the mean density 
in this region and changes the details of the gaps found 
in the region. For either centerline of the stream, the 
confidence level of the 3 to 4 gaps in the region is not 
quite 99%, but above 95% confidence. The presence of 
gaps in this region suggests that it is unlikely to be young, 
i.e. recently released from a progenitor. At the present 
time these small deviations from the smooth path of the 
stream are not very statistically significant and are not 
used in the analysis, however they are potentially a very 
powerful tool as the stream signal-to-noise is improved 
in future observations. 

4. GAP STATISTICS AND SUB-HALOS 

4.1. The Stream Age 

The mean dynamical age of the stream, that is, the 
average time since the stream stars became unbound 
from the progenitor, is a key parameter for the inter- 
pretation of the stream gaps. The ages of the stars are 
estimated from the t he photometry to be about 9 Gyr 
(jKoposov et al.ll2010f) . which sets an upper limit to the 
stream age. A simple estimate of the dynamical age of 
the stream assumes that the stars started near a common 
point (which rotates with time) but with a small spread 
in angular momentum that causes the stars to spread 
apart along nearby orbits. As the low and high angu- 
lar momentum stars drift apart, the width of the stream 
becomes entirely dependent on the dispersion of random 
velocities perpendicular to the stream orbits. The radial 
velocity dispersion is approximately fta, where a is the 
typical epicyclic size and k is the epicyclic frequency. A 
locally flat rotation curve has k — \/2Q, the local circu- 
lar velocity frequency. Taking a — 0.035 kpc at half the 
stream width, D, — 220kms~ 1 /15kpc, we find that the 
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velocity dispersion is approximately 1 k ms" 1 , well be- 
low th e 7 kms -1 upper limit measured in lKoposov et all 
(2010). For a spreading rate of 1-2 kms -1 it will take 
4.6-9.2 Gyr for stars to spread along the the 8.4 kpc long 
stream segment, if the only source of spreading was a 
range of angular momentum in the stream at some nom- 
inal time of unbinding from the progenitor. Assuming 
one end is relatively young then the mean age is 2.3-4.6 
Gyr. 

4.2. The Distribution of Gap Sizes 

The high confidence gaps in the GD-1 stream pro- 
vides an opportunity to examine the distribution of gap 
widths. The cumulative distribution of the number of 
gaps wider than some size £, N(> £), is displayed in 
Figure [5] for 99% significance gaps. The interesting fea- 
ture of the plot is that the numbers do not continue to 
rise for gaps narrower than about 1-2°. The measured 
FWHM of the s tream is 0.5°. equivalent to 70 pc. The 
ICarlbergl (|2012| ) prediction of the total number of gaps 
visible assumed that gaps as narrow as the width of the 
stream were readily visible whereas in Figure [5] the rise 
in the numbers of gaps flattens for gap sizes smaller than 
approximately 3 times the 0.5° width of the stream. 

The shallow s-curve near x = 50° of Figure [2] has gaps 
within it. Forcing the center line through this feature 
does not significantly change the distribution of gaps, al- 
though it does slightly increase the amplitude and sign- 
ficance of the features near 50°. If a stream progenitor 
disrupted at this location in the last 0.5 Gyr we would 
not expect any significant sub-halo induced gaps to have 
developed in the short tim e available. On the other hand, 
iKiipper et al.l (|2008l 12012ft have established convincingly 
that a dissolving globular cluster will have a few epicyclic 
gaps visible nearby. As a cluster dissolves its tidal ra- 
dius, which is the characteristic radius of the epicycles, 
goes to zero so that the features shri nk in proportion an d 
should become effectively invisible (JjCiippcr et al. 2010). 
Therefore a secondary conclusion is that the 50° feature, 
although not of high statistical significance in the first 
place, is unlikely to be the location of a recently dissolved 
progenitor. 

The cold-stream analysis of iCarlberd ()2012h . refined 
for galactic radii less than 30 kpc, predicts a gap creation 
rate as a function of galactocentric distance, r, in units 
of 30 kpc, and M 8 = M/10 8 M Q of, 

R u = 0.00433r°- 26 Af 8 -°- 36 kpc- 1 Gyr- 1 . (1) 

Sub-halos of mass Mg create gaps of mean length £, 

i = 9.57r°- 16 M 8 °- 31 kpc. (2) 

These relations are fits to a large set of cold stream sim- 
ulations and are only approximate. 

We eliminate the common variable between Equa- 
tions Q] and [2] to find the gap creation rate as a function 
of gap size is, 

R u = 0.060r°- 44 £~ 116 kpc^Gyr" 1 , (3) 

where the rate is the sum for all gaps larger than £, which 
is measured in kpc and the variable r is scaled to 30 kpc. 

For convenience, we convert Eq. [3] into variables ap- 
propriate to the observational quantities for the GD-1 
stream. A typical radius is approximately 15 kpc, so 
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Fig. 8. — The cumulative distribution of gap sizes at 99% con- 
fidence confidence. The lines are estimates of the expected distri- 
bution, drawn for a mean age of 4.6 Gyr (solid), 2.3 Gyr (dashed) 
and 0.5 Gyr (dotted). 
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Fig. 9. — The gap density profiles for a 70 pc wide stream en- 
countering sub-halos at 90 (solid lines) and 150 pc (dashed lines) 
from the stream centerline. The three pairs of lines of increasing 
width and depth are for masses 3x 10 s , 10 6 and 3x 10 6 Mq, respec- 
tively, for the two distances of closest approach to the centerline. 
The lines across the center are the initial density distribution and 
indicate the noise level. 

r = 0.5 and and we express gap sizes in terms of their 
measured angular sizes at 8.4 kpc, so £ = 0.140, where 
9 is in degrees. Taking the rate to be constant over the 
8.4 kpc stream segment and a mean lifetime of 4.6 Gyr 
gives, 

N{> 6) ~ 15.8t4.e6 1-1 - 16 , (4) 

which is displayed as the solid line in Figure [5J The 
relation for a 2.3 Gyr mean stream age is shown as the 
dashed line. 



GD-1 Gaps 



7 




0.1 1 10 

width [kpc] 

Fig. 10. — The relation between the estimated gap creation rate 
per unit stream length and the width of the stream. The points are 
for the streams discussed in Table [T] although only Pal 5 and GD-1 
have been analyzed with our gap-filtering methods. The hatched 
region is the theoreti cal relation for cold streams developed from 
Carlbcrg et al. (2012). The cross-hatched region assumes that the 
minimum visible gap is 3 times the stream width. For clar i ty the 
hatched regions are half the width used in Carlbcrg ct al. ( 2012) 
and underestimate the uncertainty in the theoretical relation. 

The predicted relations in Figure |8] have about the 
right slope and amplitude for gaps larger than 2°, but the 
number of gaps smaller than about 2° falls well below the 
prediction. The relationship assumes that gaps as small 
as the width of the stream are visible, whereas our gap 
analysis finds far fewer. The wa rm stream simulations of 
lYoon. Johnston fc Hoggl (|2011h chose parameters mod- 
eled on the somewhat wider Pal 5 stream, finding that 
sub-halos of mass 3 x 10 5 M Q produced gaps so small that 
they would not likely be visible at our signal-to-noise. 
Our own preliminary simulation studies of streams of fi- 
nite width find that gaps as narrow as the stream width 
are created if a sub-halo with a scale radius greater than 
the stream width passes within 1.5 times the width of 
the stream. Figure |H] shows the gap profiles in a stream 
of width 70 pc, on a circular orbit at 30 kpc radius that 
encounters sub-halos of mass 3 x 10 5 , 10 6 and 3 x 10 6 M Q 
at distances of 90 and 150 pc. These encounters are all 
for a single but reasonably representative orbit locally 
perpendicular to the stream at 1.2 times the circular ve- 
locity. At 3 x 10 5 M the gaps have a depth of ~ 30%, 
which will not be visible in the 30% noise level of the 
GD-1 stream. At 10 6 M Q the gaps are more than 60% 
deep at 150 pc and are slightly more than 30% deep at 
300 pc. Therefore these simulations support that idea 
that stream interactions effectively cut off at a mass of 
10 6 M Q and for that mass at a distance from the stream 
centerline of about 3 times the stream width. 

In young streams there is likely to be a significant range 
in angular momentum contributing to the stream width 



TABLE 1 

Observed Stream Gap Statistics 



Stream Gaps Length Width Age/2 Rgc n / n T^-u 
# kpc kpc Gyr kpc kpc -1 

Gyr" 1 



M31 


12 


200 


5 


5 


100 


6 


0.012 


Pal 5 


6 


8.1 


0.12 


3.5 


19 


22 


0.17 


EBS 


8 


4.7 


0.17 


3.5 


15 


24 


0.49 


Orphan 


2 


30 


1.0 


1.8 


30 


17 


0.037 


GD-1 


8 


8.4 


0.070 


4.6 


15 


24 


0.21 



(jEvre fc Binnevil2011[ ). Producing a gap in the stream 
requires much heavier s ub-halos, since the shear in th e 
stream blurs out the gap ([Yoon. Johnston fc Hogg|2011[) . 
The angular momentum spread at any given stream po- 
sition declines with time and will not be a problem at 
our estimated age of GD-1. 

We conclude that for the GD-1 stream the minimum 
effective sub-halo mass is 10 6 M Q which has a scale ra- 
dius of 0.11 kpc. However, depending on the orbit of 
the progenitor and the dynamics of how the stars are 
stripped away the stream material may have sufficient 
orbital shear that small gaps may be blurred out into 
invisibility. The calculation of the total number of gaps 
needs to take the distribution of gap widths into account, 
which we do with a simple empirical adjustment. 

4.3. Total Gap Rates 

Table [T] includes this paper's GD-1 measurements 
along with the star st ream data previously presented in 
ICarlberg et al.l (|2012f ) with the Pal 5 measurements. A 
simple revision to the cold stream approximation is to 
raise the minimum visible gap width to be / times the 
width of the stream, so that the minimum visible width 
w becomes fw. 

The value of / inferred from the GD-1 stream is 3 
(from a FWHM of 0.5° to an effective minimum gap size 
of about 1.5°) which reduces the coefficient of Equation[3] 
by a factor of 3.6. The value of / likely depends on the 
signal to noise in the stream and the residual angular 
momentum present along the stream, which smooths out 
gaps, but the / = 3 value should at least be represen- 
tative. The value of / = 3 corresponds to increasing 
the minimum mass sub-halo that creates a visible gap 
upwards by 3 3 22 , a factor of 35 increase, which demon- 
strates that the inferred sub-halo masses are very sensi- 
tive to both the gap measurements and the dynamical 
details of the process. 

The adjusted rate relation is displayed in Figure [TOl as 
the cross-hatched region, along with the original cold- 
stream relation as the hatched region. Note that only 
Pal 5 and GD-1 have had gaps counted with our gap- 
filtering approach. The other streams used simple count- 
ing estimates. The agreement with 99% confidence total 
counts of the Pal 5 and GD-1 streams is quite good. 

4.4. A Simple Calculation of Sub-halo numbers 

A simple model independent volumetric density for 
the objects creating the gaps, n, can be calculated from 
the basic scaling relationship for stream crossings, i?y ~ 
irnfwvrei, where v re i is the typical stream crossing ve- 
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locity and w the stream width. Taking a representative 
velocity to be the GD-1 orbit estimated circular veloc- 
ity of 220 kms" 1 (IKoposov et all 120101) . then with our 
gap creation rate of 0.21 kpc -1 Gyr -1 , gives a density of 
gap-creating objects of 0.0015 kpc" 3 for / = 3. The to- 
tal number of objects over the 30 kpc extent of the GD-1 
orbit is then about 90, where we have used the Aquarius 
density profile to compute numbers. This is quite close 
to 110 sub-halos we estimate in the same volume on the 
basis of our 10 6 M Q mass cutoff in the more elaborate cal- 
culations given above. These masses and numbers are in 
reasonable accord given the uncertainties and the strong 
dependence on the mass on the input numbers. 

5. CONCLUSIONS 

GD-1 is an exceptional star stream for the study of 
the statistics of stream gaps. The segment we examine 
has a width of 70 pc and a length of 8.4 kpc, hence a 
length to width ratio above one hundred. GD-l's sig- 
nal to noise is relatively high at about 3 per longitudinal 
bin of 0.1° (0.014 kpc) which enables a first measure- 
ment of the distribution of gaps sizes. Our study of the 
performance of our gap search procedure finds that it 
generates 44% false positives even at our 99% confidence 
level. This bias is taken into account in our gap counts, 
leading to 8 ± 3 gaps at 99% confidence. We find that 
the larger gaps are in good agreement with model predic- 
tion. However, gaps less than 3 times the stream width 
fall well below the cold stream relation. We therefore 
set the minimum visible gap size to 3 times the width 
of the stream, which is supported by some preliminary 
warm stream simulations. The revised total gap creation 
rate is in good agreement with the GD-1 data, and data 
from other streams. The total number of stream gaps 
directly indicates a population of about 100 sub-halos 
within the 30 kpc orbital extent of GD-1, which extrap- 
olates to more than 15,000 within the 433 kpc normaliz- 
ing radius of the Milky Way halo model. These numbers 
are well in excess of any known population and are in 
good agreement, within the random errors and system- 
atic error of stream age, with the distribution and total 
numbers of predicted LCDM dark matter sub-halos with 
M > 10 6 M Q . 
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